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Werner and Sornette (WS) [l| claim to show that our 
analysis of spatial distances between successive earth- 
quakes in southern California f5| does not allow to detect 
length scales associated with aftershock zones. However, 
almost all of their arguments are based on an analysis 
that is completely insensitive to finite size scaling; A 
function F obeys finite size scaling if and only if a func- 
tion / and constants a and /3 exist such that Fl{x) = 
L"- f {x / L^^) for any (large) system size L. While we stud- 
ied the probability density function (PDF) Pi(Ar) of 
spatial distances over all boxes of a given linear size L 
and its variation with L, WS do not take into account 
any variation with box size. Instead, they consider all 
events in the given catalog being in a single " box" whose 
size and shape is determined by the catalog. Since the 
scaling function / is not constrained to a particular form 
a priori, the presence or absence of finite size scaling 
cannot be established by considering a single scale L as 
claimed by WS. Thus, they cannot make any statements 
about the variation with linear size L and the existence 
of finite size scaling, which is a crucial point in our line 
of argument regarding aftershock zone scaling with main 
shock magnitude. In particular, this is true for their 
results obtained for earthquake catalogs from Northern 
California and Japan. 

In fact, the results they obtained are strikingly similar 
to results for periods of quasi-stationary seismic activity 
This suggests that the observation periods in Japan 
and Northern California that WS analyzed show rates 
of seismic activity which are rather homogeneous com- 
pared to the period we analyzed for southern California 
(1984-2000). As we emphasized in Ref. 0], the power law 
behavior of Pl ( Ar) , for earthquakes in Southern Califor- 
nia, only holds for very long observation periods where 
the rate of earthquake activity is highly heterogeneous in 
space and time. This was confirmed in a detailed study 
in Ref. The particular combination of periods and 
regions with high seismic activity and low seismic activ- 
ity in southern California is the reason for observing the 
scale-free PDF for spatial distances. Thus, considering 
only the most active period around the Landers event, 
or neglecting it, is expected to lead to strong deviations 
in the PDF as observed by WS. In particular, their re- 



sults do not in any way contradict the results we found 
and our discussion emphasizing the difference between 
homogeneous and heterogeneous rates in Q]- 

Another point raised by WS is that aftershocks occur 
at distances larger than the main shock's rupture length. 
This is well-known, yet the vast majority of what are 
typically considered aftershocks occur within distances 
which are comparable to and no more than a few times 
the rupture length (see, for example, Ref. Q and refer- 
ences therein). In particular, the largest rupture length 
in the catalog from southern California we studied is 
K, 90km. This is significantly less than the spatial ex- 
tent of the area studied (w 500km) which allows us to 
test the hypothesis of aftershock zone scaling with mag- 
nitude. As our results in Ref. ^ show, no physical length 
scale exists in the range from 20 km to « 500km. 

WS also present results for a spatially extended ver- 
sion of the ETAS model. While they claim that this 
extended model is an accurate description of aftershock 
zone scaling with main shock magnitude, a comparison 
of the Landers sequence shown in their Fig. 1 [l| and the 
ETAS model in Fig. 2 [l[ suggests otherwise. More im- 
portantly, the particular form of P{d) is speculative — it 
does not follow from the work by Kagan Q — and it is 
by no means generally accepted. Despite some indication 
for a power law with jj, k, 1.35 for short times after the 
mainshock Q, other results even directly contradict the 
form of a power-law decay for distances larger than the 
main shock rupture length 0, 0] if the activity is consid- 
ered over the long time scales relevant for our study in 
Ref. To summarize, the behavior of a model, which 
is not an accurate description of aftershock zone scal- 
ing, cannot prove that our earlier results are insensitive 
to the existence of physical length scales associated with 
aftershock zones. 
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